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The exclusive irn and pn production in hard 7*7 scattering in the forward kinematical region 
where the virtual photon is highly off-shell is studied using the 7 — » 7r _ Transition Distribution 
Amplitudes obtained in realistic models for the pion. For pn production we confirm the previous 
estimates before QCD evolution. Nevertheless, once evolution is taken into account this cross 
section grows one order of magnitude. In the case of nn production we have evaluated the cross 
section including the pion pole contribution. We observe that this contribution is responsible for 
an enhancement of two orders of magnitude with respect to the cross section evaluated without the 
0^ , pion pole term. 
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fS |. Collisions of a real photon and a highly virtual photon are an useful tool for studying fundamental aspects of 
QCD. Inside this class of processes, the exclusive meson pair production in 7*7 scattering has been analyzed in 
Ref. [l[ introducing of a new kind of distribution amplitudes, called Transition Distribution Amplitude (TDA). At 
small momentum transfer t and in the kinematical regime where the photon is highly virtual, a separation between 
the perturbative and the nonperturbative regimes is assumed to be valid. Through the factorization theorems the 
' amplitude for such reactions can be written as a convolution of a hard part Mh, the meson distribution amplitude 
<pM and a soft part, the TDA, describing the photon-pion transition, as is shown in Fig. [T] Lacking a complete 
fundamental understanding of the color dynamics, we are compelled to use models for predictions. Cross section 
5^ . estimates for the processes 

727^7r+7r- , j* L j^p + n- , (1) 

have been proposed in Ref. [2] using for the TDA a ^-independent double distributions, in a first approach, and, in a 
£^ ' second, the ^-dependent results of Ref. 0|. 

. The process 7*7 — > nn is particularly interesting because different kinematical regimes lead to different mechanisms, 
which implies a description of the process through either the pion Generalized Distribution Amplitudes or the pion- 
photon TDAs Therefore, it is of interest to deepen our understanding of the description of this process. Recently 
the pion-photon TDAs have been calculated in, respectively, the Spectral Quark Model (SQM) [5j], the Nambu - 
Jona-Lasinio model with Pauli-Villars regularization procedure (NJL) [(land a nonlocal chiral quark model 0]. A 
comparison of the results obtained in these three models is given in Ref. [8[ concluding that there is clear agreement 
between the different studies of the pion-photon TDAs, allowing us to analyze the result of a single model, e.g. the 
NJL model. Our choice is based on the fact that the NJL model Q is the most realistic model for the pion based on 
a local quantum field theory built with quarks. It gives a right description of the low energy pion physics. It has been 



o 



00 
o 



X 




FIG. 1: Factorization diagram defining the TDA for the process 7*7 — > -kM at small momentum transfer, t = (q — Pm) 2 , and 
large invariant mass, s = (p 7 + q) 2 . 
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applied to the study of pion parton distribution (PD) [l(| and pion generalized parton distribution (GPD) [ll|. In the 
chiral limit, its quark valence distribution is as simple as q(x) = 6 (x) 9(1 — x). Once evolution is taken into account, 
good agreement is reached between the calculated PDF and the experimental one [lfj. More elaborated studies of 
pion PD have been done in the Instanton Liquid Model [13] and lattice calculation based models [l3[ using nonlocal 
Lagrangians 14j, which confirm that the result obtained in the NJL model for the PD is a g ood approximation. The 
QCD evolution of the pion GPD calculated in the NJL model has been also studied in [la ]. 

Before defining the TDAs, we introduce the light-cone coordinates = (v° ± v 3 ) j\pl and the transverse com- 
ponents v 1 - = (u 1 ,^ 2 ) for any four-vector v^. We define P = (p n + Pj) /2 and we introduce the light-front vec- 
tors = P+ (1,0,0,1) /y/2 and tt/ 1 = (1,0,0,-1) / (v / 2P+) . The momentum transfer is A = p„ - p 7 , therefore 
P 2 = m 2 r /2 — t/A and t = A 2 . The skcwness variable describes the loss of plus momentum of the incident photon, i.e. 



£ = ( Pl — Ptt) + /2P + , and its value ranges between — 1 < £ < 
axial TDA are defined by 



-tj (2m 2 — t). With these conventions, the vector and 
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where the pion decay constant is / w = 92.4. MeV, e (£) is equal to 1 for £ > and to —1 for £ < and (/> n (x) is the 
pion DA. Here we have modified the definition given in Refs. [3, H| m order to introduce the pion pole contribution 
0, Q in Eq. ([3]). This equation deserves some comments. The pion pole term in ([3]) describes a point-like pion 
propagator multiplied by the distribution amplitude (DA) of an on-shell pion. It contributes to the axial current 
through a different momentum structure and must be subtracted in order to obtain de axial TDA. We emphasize that 
it is a model independent definition, because we have define the numerator of the pion pole term as the residue at the 
pole t = m 2 . With this definition, all the structure dependence related to the outgoing 7r =b is included in A(x, £,t). 
Moreover, the pion pole contribution can be estimated in a phenomenological way, as we will see later on. With these 
definitions we recover the sum rules 



dx D(x,£,t) = ^lll- F D (t) 



D = V,A 



(4) 



with the standard definitions for the form factors Fv,a appearing in the tt ± ->■ ^v~i decay [l6[. Notice that the 
on-shell pion DA obeys the normalization condition J dx<p 7r (x) = 1. 

The 7*7 — > M + n~ process, with M + = p\ or 7r + , is a subprocess of the e ip e )+l (Py) — * e {p' e ) + M + (pm) + k~ (Ptt) 
process. We follow all the definitions of the kinematics given in the section III. A and Fig. 3 of Ref. 0], with the 
exception that our vector is twice the nP" vector used in Q (i.e. n.p = 1 with our definitions). In particular, for 
massless pions, we have 



Q 



(Pe -PeY 



(Pe+PfY 



(5) 



Py={l + 0P > P* = (l~flP+ 2(l-fl Tt + ^ ± ' 

q=-2&+^n , (6) 

where A T = (0,^,0) and therefore A 2 , = -A ±2 . Notice that A ±2 = (-t)(l - £)/(l + £), with t < 0. The 
longitudinal polarization of the incoming virtual photon is Sl — (2£,P/Q + Q n /{^€))- The real photon polarization is 
defined by the condition e~ = together with the gauge condition e + = 0. 
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The differential cross section is given by 1 0] 



with 



v ~ ~q ] dx I dz — 4> p {z) 
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(7) 



(8) 



(9) 



where z is the light-cone momentum fraction carried by the quark entering the meson M + and f p = 0.216 GeV and 
Q q is the electric charge of the quark q. The last term on the r.h.s. of Eq. ([9]) is the pion pole contribution to the 
amplitude coming from the second term of Eq. ([3]). 

We proceed now to evaluate these integrals. The meson DA is chosen to be the usual asymptotic normalized meson 
DA, i.e. 4>m(z) — 6z(l — z), which cancels the z-dependence of the hard amplitude. For the nonpcrturbative part 
of the process we use the TDAs evaluated in the NJL model. This approach is based on the determination of the 
pion as a bound state through the Bethe-Salpeter equation. This guarantees that all invariances of the problem are 
preserved. As a consequence, the obtained TDAs explicitly verify the sum rules, the polynomiality condition, the 
isospin relations and have the correct support in x [6j . 

In Ref. @ the ir + — > 7 TDAs were calculated, which are connected to the 7 — > ir~ TDAs through CPT symmetry 



D~ 



D = V 1 A 



(10) 



In Figs. [2] and [3] are shown the vector and axial 7 — * ir~ TDAs for t = —0.5 for different £ values. From Eq. ([7|) 
it can be observed that £ > Q 2 / (2s e7 — Q 2 )- In other words, there is a (positive) lower limit on the value of £. It is 
indeed a particularly interesting restriction because the value of £ defines the shape of the TDAs. In particular, the 
shape of the axial TDA radically changes according to the sign of the skewness variable; A(x,£,t) has, at x = ±£, 
its maximum values for £ > (see Fig. [3]) while it has its minimum values for £ < [8] . However the vector TDA 
has its maximum and minimum values at x = ±£ (see Fig. [2]) independently of the sign of the skewness variable. On 
the other hand, the magnitude of the distributions is controlled by the i-dependence. This can be easily understood 
because the TDAs, that must satisfy the sum rules Eq. ([4|), are expected to decrease at least as t~ l . 

In order to numerically estimate the cross sections, we need to fix the strong coupling constant. In Ref. [li| is 
mentioned that a large value of a s (a s ~ 1) is indicated in the case of an asymptotic DA. Using this value for a s 
we have evaluated the cross section for p production. In Fig. |4] we plot this cross section as a function of £. As we 
observe, the cross section is largely dominated by the imaginary part of the integral of Eq. ([5]) . The t-dependence 
of the cross section comes from both the overall (— t) factor present in Eq. ([7]) and the t dependence of the Vector 
TDA. Therefore, we expect a decreasing as i" 1 for large t values. Comparing with the previous results in Ref. we 
observe that our predictions are higher by a factor 2 or 3. 

The 7r production is described through Eq. (|9|). Here we have a pion pole term which becomes proportional to the 
electromagnetic pion form factor (FF) 



dx 



Qa 



A 1 ' 



3 Q 2 F„ (Q 2 



4tt t 



(11) 



If we use the asymptotic form of the pion DA, i.e. W (z) with z = (x + £) /2£, in Eq. 
pion FF 



we obtain the Brodsky-Lepage 



Q 2 F„ (Q 2 ) = 16tt a s f 2 



(12) 



1 A factor 1/4 is missing in Eq. (23) of Ref. 0]. This typo does not affect to the numerical results reported in that Reference [17| . 
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FIG. 2: The functions V 7 ^ (x,£,t) and for £ = 0.25, 0.5, 0.75 and for t = -0.5 GeV 2 . In each figure, the solid line 
corresponds to the NJL model prediction and the dashed line to its LO evolution. 
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FIG. 3: The functions A 1 ^ 7 ' {x,^,i) for different values of the skewness variable £ and for t = —0.5 GeV 2 . 



The cross section for pion production at Q 2 = 4 GeV 2 as a function of £ is given in Fig. [5]4eft). This cross section is 
dominated by the pion pole contribution which is determined by the Brodsky-Lepage pion FF. Alternatively, if the 
pion FF is experimentally known, we can infer phenomenologically this contribution. And hence the axial TDA could 
be extracted from the interference term. From Ref. [20j we know that the pion FF at Q 2 = 2.45 GeV 2 is 0.167 ± 0.010. 
In Fig. fright) we depict, for each contribution, the prediction using the interval defined by the experimental value 
of the FF (filled areas), including also the theoretical prediction using the Brodsky-Lepage pion FF (lines). 

The ^-dependence of the cross section for pion production includes a strong dependence on t coming from the 
pion pole. Neglecting the pion mass in ([9]), we observe that the pion pole contribution to X v is proportional to t . 
Therefore, the cross section grows as i _1 for small t values. For large t values we expect, as in the p production case, 
a decreasing as t~ x . 

We have also studied the effect of the QCD evolution on our estimates, using for this purpose the code of Freund 
and McDermott 2 if . One needs to fix the value of Qo for which the quark distributions obtained in the NJL model 
are considered to be a good approximation of the QCD quark distributions. Knowing that the momentum fraction of 
each valence quark at Q = 2 GeV is 0.23 [22], we fix the initial point of the evolution in such a way that the evolution 
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FIG. 4: e-y — > e'p^-K^ differential cross section plotted as a function of £ for Q 2 = 4GeV 2 , s e7 = 40GeV 2 , t = -0.5 GeV 2 . In 
the first layer, the dotted (dashed) line is the contribution to the cross section coming from the imaginary (real) part of the 
integral given in Eq. In the second layer we give the cross sections before and after evolution. 




FIG. 5: cy — > e'-R + -K differential cross section plotted as a function of £ for s el = 40 GeV 2 , t = -0.5 GeV 2 , Q 2 = 4 GeV 2 (left) 
and Q 2 = 2.45 GeV 2 (right). The dashed (dashed-dotted) [dotted] line is the contribution to the cross section coming from the 
pion fom factor (axial TDA) [interference term between the pion FF and A]. The pion pole contribution is calculated using the 
Brodsky-Lepage pion FF. The filled areas in the right layer correspond to the same contributions but with the experimental 
value for the pion FF F n = 0.167 ± 0.010 



of the second moment of the pion parton distribution reproduces this result. This condition is fulfilled at a rather 
low value, i.e. Qo — 0.29 GeV, when the LO evolution is used. Going to the NLO changes this value to Qq = 0.43 
GeV. However, in the latter case, the resulting PD is basically unaffected by the change in Q. The effect of the NLO 
evolution is compensated in the LO evolution going to a lower value of Qo, a result that has already been noticed in 



proton parton distributions 23]. In order to illustrate the latter statement, we have depicted both the evolved pion 



PD for the LO and the NLO in Fig. [6] Turning our attention to the vector TDA evolved at LO (Fig. [2]), we observe 
that the value of V (x, £, t) at x = ±£ grows for small values of £ and decreases for large £ values in comparison with 
the TDA at the scale of the model. This implies that the cross section for p production, which is largely dominated by 
the imaginary part, will grow appreciably in the small £ region. In Fig. |4]we compare the cross section after evolution, 
calculated only through contribution of the imaginary part of X p , with the one obtained before evolution. We observe 
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FIG. 6: Pion parton distribution. The solid line corresponds to the LO evolution of the NJL prediction and the dashed line to 
the NLO. Experimental data are from (24j|. 



that this cross section is multiplied by a factor about 5 in the £ ~ 0.2 region. 2 In the case of the axial TDA, the cross 
section is dominated by the pion FF contribution, therefore the effect of the evolution is expected to be small. 

In many papers present in the literature, the r.h.s. of Eq. ([3]) contains only the A term. In a general case the A 
term and the pion pole terms have different tensor structures, but we can fix the gauge convention in such a way that 
these two structures coincide, as we have done in this paper. In that case we can change our definition of the axial 
TDA 

/ ^ ixP+z ~ <^M*(-f) 7+ '» T±q (D hM) 

with 

(z, £, t) = A^ (a, £, t) + e{£) fa (^f) . (14) 

mi-t V 2 W 

The latter expression shows that the pion pole contribution to the axial TDA is closely related to the D-term of 
the Generalized Parton Distributions [25|. Nevertheless, it must be realized that in this case it gives an explicit 
contribution to the sum rule 

/ dx A^ (x, £ t) = F D (t) + iPy - P*) -n . (15) 

In this paper we have looked at the expected cross section for 7r-7r and n-p production in exclusive 7*7 scattering 
in the forward kinematical region using realistic models for the description of the pion. First we confirm the previous 
estimates for p production, even if our results for the cross section are a factor 2 larger than the one obtained in 
Ref. [2|. Second, in comparison with this previous evaluation of the cross sections, we have improved in considering 
the effect of evolution on the vector current. In doing so, an additional factor 5 appears in the small £ region, leading 
to a cross section for p production of one order of magnitude larger than the previous calculation. We have also 
improved in including the pion pole term in the tensor decomposition of the axial current. Then an even larger 
enhancement factor, of about 60 in this case, is found in the cross section for pion production. The interference term 
becomes a factor 15 larger than the pure TDA contribution, making the axial TDA more accessible experimentally. 



A 1 ' 



V2U 



(13) 



2 The effect of QCD evolution, if calculated also through the contribution of the real part of T p , could lead to a change in the cross section 
of about 15%, what is within the model's uncertainties. 
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The interest on TDAs is actually extended to other transitions such as 7*iV — > N'it, j*N — > iV'7, iVTV — > 7*7, which 
are related to AT-7T and iV-7 transition distribution amplitudes 26|, 27, 28j |. 

We are thankful to J. P. Lansberg for useful discussions. This work has been supported by the Sixth Framework 
Program of the European Commission under the Contract No. 506078 (13 Hadron Physics); by the MEC (Spain) 
under the Contract FPA 2007-65748-C02-01 and the grant AP2005-5331 and by EU FEDER. 



Electronic address: Aurore.Courtoy@uv.es 
Electronic address: Santiago.Noguera@uv.es 

B. Pire and L. Szymanowski, Phys. Rev. D 71 (2005) 111501 [arXiv:hep-ph/0411387|. 

J. P. Lansberg, B. Pire and L. Szymanowski, P hys. Rev. D 73 (2006) 074014 [arXiv:nep-ph/0602195l . 

B. C. Tiburzi, Phys. Rev. D 72 (2005) 094001 [arXiv:hep-ph/0508112j. 

I. V. Anikin, I. O. Cherednikov, N. G. Stefanis and O. V. Teryaev, LarXiv:0806.4551 [hep-ph]. 
W. Broniowski and E. R. Arriola, Phys. Lett. B 649 (2007) 4 9 [arXiv:hep -ph/0701243 . 

A. Cou rtoy and S. N oguera, Phys. Rev. D 76 (2007) 094026 |arXiv:0707 33661 [hep-ph]]. Prog. Part. Nucl. Phys. 61, 170 

(2008) [arXiv:0803.3524l [hep-ph]]. 

P. Kotko and M. Praszalowi cz, larXiv:0 803.2847 [hep-ph]. 
A. Courtoy and S. Noguera. larXiv:0 804.4337 [hep-ph]. 

Y. Nambu and G. Jona-Lasinio, Phys. Rev. 122 (1961) 345. . S. P. Klevansky, Rev. Mod. Phys. 64 (1992) 649. 

R. M. Davidson and E. Ruiz Arriola, Phys. Lett. B 348 (1995) 163. Acta Phys. Polon. B 33 (2002) 1791 

[arXiv:hep ph/0110291] . 

S. Noguera, L. Theussl and V. Vento, Eur. Phys. J. A 20 (2004) 483 [arXiv:nucl-th/0211036j . 

I. V. Anikin, A. E. Dorokhov and L. Tomio, Phys. Part. Nucl. 31, 509 (2000) [Fiz. Elem. Chast. Atom. Yadra 31, 1023 
(2000)]. 



S. Noguera and V. Vento, Eur. Phys. J. A 28 (2006) 227 arXiv:hep-ph/0505102 
S. Noguera, Int. J. Mod. Phys. E 16 (2007) 97 |arXiv:hep-ph/0502171| .~ 
W. Broniowski, E. R. Arriola and K. Golec-Biernat, Phys. Rev. D 77 (2008) 034023. [arXiv:0712.ToT2l [hep-ph]] 
C. Amsler et al. [Particle Data Group], Review of particle physics, Phys. Lett. B667 (2008) 1. 
J. P. Lansberg, Private Communication. 



V. M. Braun, A. Lenz, G. Peters and A. V. Radyushkin, Phys. Rev. D 73 (2006) 034020 [arXiv :hep-ph/0510237]. 
G. P. Lepage and S. J. Brodsky, Phys. Lett. B 87 (1979) 359; SLAC-PUB-2294; A. V. Efremov and A. V. Radyushkin, 
Phys. Lett. B 94 (1980) 245. 

T. Horn et al [Jefferson Lab F(pi)-2 Collaboration], Phys. Rev. Lett. 97, 192001 (2006) [arXiv:nucl-ex/0607005] . 
A. Freund and M. F . McDermott, Ph ys. Rev. D 65, 056012 (2002) [Erratum -ibid. D 66, 079903 (2002)]. 
arXiv:hep-ph/0106115 . .See the web page http://durpdg.dur.ac.uk/hepdata/dvcs.html 
P. J. Sutton, A. D. Martin, R. G. Roberts and W. J. Stirling, Phys. Rev. D 45 (1992) 2349. 
M. Traini, A. Mair, A. Zambarda and V. Vento, Nucl. Phys. A 614 (1997) 472. 
J. S. Conway et al., Phys. Rev. D 39 (1989) 92. 

M. V. Polyakov and C. Weiss, Phys. Rev. D 60 (1999) 114017 [arXiv:hep-ph/99024511. 

J. P. Lansberg, B. Pire and L. Szymanowski, Nucl. Phys. A 782 (2007) 16 arXiv :hep-ph/0607130| . 

J. P. Lansberg , B. Pir e and L. Szymanowski, Phys. Rev. D 75, 074004 (2007) [Erratum-ibid. D 77, 019902 (2008)] 

[arXiv:hep-ph/0701125] . 

[28] M. Pincetti, B. Pasquini and S. Boffi. larXiv :0807.4861 [hep-ph]. 



